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FOREWORD 

This  report  was  prepared  In  the  Mechanics  and  Surface  Interactions 
Branch  (AFWAL/MLBM),  Nonmetalllc  Materials  Division,  Materials  Laboratory, 
Air  Force  Wright  Aeronautical  Laboratories,  Wrlght-Patterson  Air  Force 
Base,  Ohio.  The  work  was  performed  under  the  support  of  Project  Number 
2307,  "Nonmetalllc  Structural  Materials",  Task  Number  2307P2,  "Composite 
Materials  and  Mechanics  Technology". 

In  this  report,  an  automated  composite  laminate  sizing  technique  Is 
presented,  which  optimizes  for  minimum  weight.  The  technique  can  be 
coded  for  a  microcomputer,  and  a  listing  In  BASIC  is  given  for  an  Epson 
HX-20  microcomputer.  The  program  is  interactive  and  easy  to  use.  The 
optimization  Is  for  point  stress  under  multiple  loads. 

The  program  Is  available  on  an  audio  cassette,  and  can  be  obtained 
by  sending  a  blank  tape  to  AFWAL/MLBM,  Wrlght-Patterson  AFB,  Ohio  45433 
and  referencing  this  report. 
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SECTION  I 

PROGRAM  DESCRIPTION 

CLASS  (Composite  Laminate  Automated  Sizing  for  strength)  Is  an 
Interactive  optimization  program  designed  to  run  on  a  small  microcomputer. 
The  listing  presented  here  Is  for  an  Epson  HX-20  portable  computer  with 
16K  of  memory.  The  version  of  Basic  Is  standard  enough  that  translations 
to  other  microcomputers  would  be  easy. 

The  program  will  find  a  minimum  thickness  laminate  which  will  not 
fail  under  any  of  the  load  conditions  entered.  Ply  orientations  are 
chosen  by  the  user.  The  program's  capability  in  handling  multiple, 
independent  loads  could  be  useful  for  loads  which  change  with  time  or 
position  on  a  constant  thickness  plate,  or  for  situations  where  there  is 
uncertainty  in  calculating  the  loads.  As  the  program  is  currently 
dimensioned,  four  Independent  load  combinations  and  6-ply  orientations 
can  be  entered. 

Only  point  stresses  are  considered,  thus  the  program  optimizes  the 
laminate  only  at  one  point  in  the  structure.  Furthermore,  the  program 
assuems  lrplane  loads  only  and  no  out-of-plane  deflections.  This  implies 
a  symmetric  laminate,  but  stacking  sequence  is  not  a  factor  in  the 
program.  The  layer  thicknesses  generated  by  the  program  are  the  total 
and  must  be  divided  by  2  to  get  the  halves  of  a  symmetric  laminate. 

No  knowledge  of  optimization  techniques  is  needed  to  run  the  program 
and  very  little  knowledge  of  laminate  plate  theory  is  needed.  In 
addition,  material  properties  for  five  common  advanced  composites  are 
stored  in  the  program,  or  the  computer  can  ask  for  new  properties  through 
prompts. 
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I 

\ 

\  SECTION  II 

| 

|  GENERAL  INSTRUCTIONS 

I 

The  program  can  be  entered  by  hand  from  the  listing  given,  or  by  an 
audio  cassette  tape  available  from  AFWAL/MLBM.  If  entered  by  hand, 
material  properties  can  be  loaded  by  running  the  program  and  using  the 
"NEW  material  option.  To  load  the  tape,  simply  hook-up  to  cassette  as 
shown  In  the  Epson  Instruction  manual,  and  enter  LOAD  "CLASS",  R.  The 
"R"  is  needed  to  automatically  start  a  part  of  the  program  which  In  turn 
loads  the  material  properties  from  the  tape.  Material  properties  are 
stored  In  the  Epson's  "Ram  File"  feature  which  allows  the  machine  to  be 
turned  off  without  losing  data. 

Running  the  program  on  another  computer  should  be  possible  If  13K 
of  memory  Is  available.  The  commands  most  likely  to  need  changing  are 
PUT*  and  GET*  used  for  the  Epson's  Ram  File  operations.  Equivalent  disk 
commands  or  data  statements  can  be  used  Instead. 

When  running,  the  program  prompts  the  user  for  Input  Information. 

An  example  dialogue  between  computer  and  user  Is  given  below  along  with 
the  printed  output  generated. 
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SECTION  III 

COMPUTER/USER  DIALOGUE 


LCD  Display 


Keyboard  Response 
(comments  in  parenthesis) 


Press  any  key  when  desired 
Material  appears 

T300/5208 

B  ( 4 ) /5505 

AS/3501 

Scotchply  1002 

Kevlar  49/Epoxy 

Aluminum 

New 

REVIEW  OR  NEW  DATA  (R/N)  ? 

WHICH  MATERIAL  WILL  YOU 
REPLACE  (0-5)  ? 

EX(GPa)  =  ? 

EY(GPa)  =  ? 

VX  =  ? 

ES(GPa)  =  ? 

X (MPa)  =  ? 

X’ (MPa)  =  ? 

Y (MPa)  =  ? 

Y' (MPa)  =  ? 

S (MPa)  =  ? 

THICKNESS  (m. )  =  ? 

NAME  (15  CHR  MAX)  ? 

ADDITIONAL  CHANGES  ? 


RUN  RETURN  (unless  otherwise 
noted,  "Return"  key  pressed 
after  each  keyboard  entry) 


(random  key  pressed  when 
"New"  appears  on  screen) 

N 

5 

(materials  numbered  in  same 
order  as  listed  T300/5208=0) 

185 
6.76 

.2 

5.86 

680 

690 

(primed  constants  imply 
compressive  properties) 

16 

186 
72 

125E-6  (ply  thickness) 
HMS/3002M 

N 
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LCD  Display 


HOW  MANY  PLY  GROUPS  ? 
ENTER  PLY  GROUP 
ORIENTATIONS 
PLY  1  =  ? 

PLY  2  =  ? 

PLY  3  =  ? 

PLY  4  =  ? 

ENTER  NUMBER  OF  INDEPEN¬ 
DENT  LOADING  CONDITIONS  ? 


LOAD  1 

in  MPa 

Nl 

=  ? 

N2 

=  ? 

N6 

=  7 

LOAD  2 

in  MPa 

Nl 

=  7 

N2 

—  ? 

N6 

=  7 

WORKING  ITERATION  1 

TOTAL  THICKNESS  = 

1.71342  E-02  m. 

137.07  PLIES 
HIT  ANY  KEY  TO  CONTINUE 

Press  Y  if  printout  of  displayed 
result  is  desired.  Press  N  if 
not 


Keyboard  Response 
(comments  in  parenthesis) 

(Materials  list  begins  again, 
this  time  with  the  new  material 
replacing  aluminum,  when  it 
appears  a  key  is  pressed) 

4 


0 

90 

45 

-45 


2 


3 

2 

.5 


1 

4 

0 

(after  4  iterations  and  about 
7  minutes  the  computer  beeps 
that  the  solution  has  been 
found.  This  example  ran  for 
an  unusually  long  time.  Most 
problems  will  run  in  less  time) 

(press  any  key,  no  return) 
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Material  Properties 
HMS/3902M 
EN=  185  GPa 
EV=  6.  76  GPa 
ES=  5.86  GPa 
UX=  .2 
X=  680  MPa 
X’=  698  MPa 
V=  16  MPa 
7’=  186  MPa 
S=  72  MPa 

Ply  Thickness  .080125  m 


LOADING  1 
N  1=  3  MN'  fri 
N  2=  2  MN-'m 
N  6=  .5  MN'  Ui 
LOADING  2 
N  1=  1  MN/m 
N  2=  4  MN-m 
M  6=  9  MN/m 


Norm.  I A|  in  GPa. 


i  74.7591 

6.  51 1  i 

5.  549 

!  6. 5111 

106. 9691 

5.  549 

!  5.  549  i 

5.  5491 

11.  017 

Compliance  C normalized) 
in  1/TPa. 

i  13.  922  i 

-8. 4971 

-6.  762 

i  -8.4971 

9.  6171 

-4.  594 

!  -6.7621 
*  -  * 

-4.  594 i 

-  -  -  ‘ 

96.  490 

**Ply  Patio  Graph** 
ANGLE 


T  ota  1 

thickness= 

+0MMMI 

.  0171E+00  Hi. 

137.  i. 

)7  Plies 

ANGLE 

RATIO 

#PL I ES 

+45  m 

0 

.  3476 

47.  65 

-45 

1  1  1  1  1  1 

98 

.  5281 

72.  38 

45 

.  1243 

17.  04 

t 

-45 

0 

0 

.  10 

PLV  RATIO 

STRENGTH  RATIOS 
1  ULTIMATE  STRAIN 
>1  IS  SAFE 


LOADING 

1 

PL  ' 

0 

1.  1622 

90 

1 

45 

1.  1616 

LOADING 

2 

PLV 

0 

1.  0054 

98 

1.  2115 

45 

1. 0122 

LAMINATE  STRAINS 
LOADING  1 
el=+2.  1S2E-83 
e2=+0.  902E-83 
e6=+l.  O96E-03 
LOhDING  2 
el=+0.  697E-03 
e2=+2.  216E-83 
e6=-l.  467E-03 


i 

i 


! 


i 

j 


Figure  1.  Output  Produced  from  Example  Input  Dialogue 
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LCS  Display 


Keyboard  Response 
(comments  in  parenthesis) 


PLY  PROPERTIES 

LOADS 

TOTAL  THICKNESS  & 
PLY  RATIOS 

STRENGTH 

LAMINATE  STRAINS 
STIFFNESS  MATRIX 
COMPLIANCE  MATRIX 
PLY  RATIO  GRAPH 

FINISHED 

HIT  ANY  KEY  TO 

CONTINUE 


Y  (return  key  not 
used  for  these  responses) 

Y 


Y 

Y 

Y 

Y 

Y 

Y 

(after  entire  list  of  print 
out  options  is  presented, 
computer  produces  the  print 
out  shown  on  next  page) 

(pressing  a  key  restarts 
program.  Press  "BREAK” 
key  to  exit) . 
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SECTION  IV 
METHOD 

The  goal  is  to  minimize  the  total  thickness  of  a  composite  laminate 
subject  to  failure  constraints  under  static  loads.  Specifically, 

m 

2  h.  =  min  where  m= number  of  ply  groups 
1=1  ' 

subject  to  h^  ^  0 


(ej  (L)  (L)  (e  )  (L) 

and  G..  e.  e.  +  G.  e.  -  1  0  where  h.  is  the  total 

J  K  J  K  J  J  1 

thickness  of  all  the  plies  at  the  kth  orientation  (which  will  be  referred 
to  as  a  "layer"  in  this  report).  The  failure  criteria  is  a  first  ply 
failure  based  on  the  Tsai-Wu  tensor  criteria  in  strain  space.  The  G's 
are  transformed  to  the  laminate  axis  from  the  ith  layer's  orientation. 

The  strains  are  associated  with  the  Lth  loading  combination.  This 
distinction  is  made  since  more  than  one  independent  loading  may  be  con¬ 
sidered.  For  the  definition  of  the  G's  in  terms  of  experimental  strength 
data,  see  Reference  1. 

Stacking  sequence  is  not  included  in  this  formulation,  and  the 
laminate  is  assumed  not  to  bend  or  warp.  Therefore,  strains  and  loads 
are  related  by 

7T=.ai7 

The  optimization  method  applied  is  a  modification  of  the  method  of 
feasible  directions  (Reference  2).  The  method  can  be  demonstrated 
graphically  with  two-dimensions,  i.e.  two  layers.  In  Figure  2  the  two 
equalities 


(0)  (0) 

V.VGI  v° 

a«.°>  <901,  _,.o 

IJ  I  I  I  I 


7 


>*  y  •» 


- v-1  '  1  *  * 

ti  m >iifrA — 
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have  been  ploted  as  functions  of  h^  and  h^-90^  for  the  single  loading 
condition  shown.  Any  point  above  and  to  the  right  of  these  two  curves 
Is  feasible,  that  is,  failure  will  not  occur.  Points  to  the  left  and 
below  the  curves  are  infeasible.  Because  our  objective  function  (the 
sum  of  the  layer  thicknesses)  is  linear,  the  optimum  point  will  lie  on 
one  of  these  curves  or  the  Intersection  of  multiple  curves. 


2  9  4  5 


Figure  2.  Constraint  Surfaces  and  Optimization  Trajectory  for  0/90 
Under  Biaxial  Load 
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The  program  starts  by  finding  an  initial  feasible  point  (A)  which 
lies  on  a  constraint  curve  farthest  from  the  origin  on  the  line 
h^^  -  The  distance  from  the  origin  is  calculated  using  a  strain 

ratio  method.  Along  any  vector  which  passes  through  the  origin,  new  ply 
group  thickness  coordinates  are  scaled  by 

h,  =  h®  •  S/$° 

where  S  is  a  scalar  distance,  h?  are  the  coordinates  of  the  current  point 
and 

s °.[2  (h®)»]w*  . 

li»i  1  j 

Along  this  vector,  strain  can  be  found  using 


where  e°  is  a  component  of  laminate  strain  evaluated  at  SQ.  Substituting 
into  the  failure  criteria  we  have 


jQi)  o  o  o2 

6i«  ♦)  «>  s 

r-02 


G^’.»  S° 


S02 


-  I«0 


To  ensure  the  calculated  point  lies  slightly  in  the  feasible  region 
despite  any  numerical  error,  the  program  sets  this  function  equal  to  the 
negative  of  a  small  number  (e^)  rather  than  zero.  Solving  this  equation 
for  positive  S  we  have 


-B+ZCB*-  4AC) 
2A 


9 


r  M  •  %*- 
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where 


A»l-e 


*  (8,) 

j*i  j  J 


.  A  O 

B*  Z-6  «,  S 


c»£  £-«,‘?l’v»sM 

jk  J  " 


J=l  k=l 


If  S°  lies  in  the  feasible  region  we  solve  the  above  equation  for  each 
layer  and  each  load  combination,  and  then  take  the  smallest  resulting  S 
as  the  one  that  defines  the  boundary  of  the  feasible  region. 


The  next  step  in  the  optimization  procedure  is  to  establish  a 
direction  vector  which  will  point  away  from  the  constraint,  point  A  lies 
on  and  is  parallel  to  the  plane  defined  by  Eh^  =  constant.  In  Figure  2, 
this  direction  Is  shown  as  Z.  Finding  Z  first  requires  calculation  of 
the  gradient  of  the  active  constraint  evaluated  at  A.  Let 


:P,L* 


6 


(Bp)  (l)  (U  (Bp)  (L) 
')»*) r  - 


+  6. 

i 


«. 

I 


-  I 


where  k  and  N  correspond  to  the  layer  and  load  combination  of  the  active 
constraint.  A  constraint  is  considered  active  if 

cP,L*-», 

where  e2  is  a  small  number.  Note,  that  more  than  one  constraint  may  be 
active.  The  gradient  is  then  given  by 


-»  m 
VCB 


CO 


(9p)  dt 


J* 


)  +0 


(0p)  d* 


<L) 


A 


where  h^  Is  a  unit  vector.  To  find  the  partlals  of  strain,  we  start 
with  the  basic  equation 


N*  IAI  « 


10 


* 
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and 


o*/-|a|7  +  a/-7 

dh,  dh, 

£~~'  *  “I  a  1 1  IaI  7 

dh,  dh, 


dh( 


lAl 


(01) 

(01) 

(01) 

Qu 

«I2 

Q!3 

(0i) 

(00 

(0i) 

°2I 

°22 

(0i) 

(01) 

(00 

°I3 

°23 

°33 

rJ0l . 
*  [Q,  1 


The  gradient  vector  is  normalized  to  unit  length.  If  more  than  one 
constraint  is  active,  the  normalized  gradients  are  sumned  together  and 
the  sum  is  then  normalized  to  one.  The  negative  of  the  gradient  will 
point  away  from  the  constraint,  into  the  feasible  region.  This  vector 
is  now  projected  onto  the  plane  defined  by  the  unit  normal  n,  where 


a  l  2  a 

Aw  1 


The  projection  can  be  made  with  a  double  cross  product 


1 


A 

*n 


x  (- Vc  i  n ) 


With  a  vector  identify,  this  can  be  rewritten  as 

Z*(VCn)n-Vc 

Finally,  Z  is  also  normalized  to  unit  length. 

Along  Z,  another  constraint  will  eventually  be  reached  (point  B  in 
Figure  2).  The  point  is  found  iteratively  by  a  bisection  technique. 
Since  the  bisection  method  is  very  time  consuming,  the  constraint  line 
is  only  found  within  a  relatively  large  error  band.  What  we  are  really 
interested  in  is  a  point  approximately  midway  between  A  and  B,  which  is 
C  in  the  figure.  From  point  C,  the  strain  ratio  technique  is  used  to 
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analytically  calculate  0.  Starting  at  D,  the  entire  procedure  repeats. 
The  program  terminates  when  the  distance  AB  or  CD  is  small  (say  1/10  a 
ply  thickness)  or  the  magnitude  of  Z  before  normalization  is  very  small 
(implying  h  and  vc  are  almost  parallel). 

In  some  cases,  h^  >_  0  constraint  may  be  reached.  When  this  happens, 
that  orientation  is  completely  dropped  from  further  calculations.  Thus, 
the  constraints  associated  with  a  zero  thickness  layer  cannot  effect  the 
results.  Once  an  orientation  reaches  zero  thickness.  It  is  never 
reinstated  In  later  iterations. 

Figure  1  shows  a  case  where  the  program  reaches  the  Intersection  of 
two  constraints.  However,  simultaneous  failure  should  not  be  considered 
a  criteria  for  optimization.  Figure  3  shows  a  case  where  only  one  layer 
approaches  failure.  The  constraint  line  for  the  +45°  layer  Is 
completely  in  the  infeasible  region.  The  line  h^5^  +  h^--*5-*  *  constant 
has  been  included  to  show  that  point  D  is  the  minimum  thickness. 


Figure  3.  Constraint  Surfaces  and  Optimization  Trajectory  for  +  45 
Laminate  Under  Pure  Shear 
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Comments 


10  •**  MAIN  CLASS** 

20  CLEAR  75.330 
30  WIDTH  21.5 
40  DEFFIL  55.0 
50  DEFIHT  I-R.'DEFDBL  F 
60  DIM  A<3» 3> > B(b> 9) >C(6 
. 6>.D<3. 3>.G<3.3)» XN<4»  3 
>. AI<3.3)»Q<3»3)»HT6>.R< 
3>,S<3>»T<6>»IK5>.U<7},X 
(6).V<3>.2<6>.E<4,3> 

70  DIM  C’-.<  10.2) 

80  DEF  FND£G<X>=X*57.  285 
73 

83  DEF  FNRAD<X>=X/57.  285 
78 

105  RESTORE 

110  READ  IMAX.E2.E5.E6 

120  ITER=1 

130  GOSUB  2540 

140  PRINT“CORRECTIONS":I 

NPUT 

150  IF  A**"V"  THEN  130 
155  CLS!  PR I NT" WORKING": 
PR I NT " I TER AT  I ON I TER 
170  GOSUB  2990 
180  GOSUB  2330 
IS©  GOSUB  2180 
200  GOSUB  1680 
205  CLS!  PR I NT "WORKING": 
PRINT" ITERATION"} ITER 
210  IF  FI=“FAIL"  THEN  33 
00 

220  GOSUB  1370 

230  ITEP=ITER+1 

240  IF  Fi=“FAIL"  OR  ITER 

>IMAX  THEN  3300 

250  GOTO  200 


20-40  commands  to  configure  the 
machine 


50  Implicit  integer  and  double 
precision 

80-90  convert  radians  to  degrees 
and  degrees  to  radians 

130  -  Gosub  input 

170  -  Gosub  invariants 

180  -  Gosub  transformations 

190  -  Gosub  initial  feasible  pt. 

200  -  Gosub  direction 

220  -  Gosub  new  thickness 


260  ’**  CONSTRAINT  TEST* 
* 

270  G*="PASS"«  NC=0 

230  FOR  P=1  TO  MPLV 

280  IF  H<P>=0  THEN  445 

300  II=P  : GOSUB  1230 

310  FOR  N=1  TO  NL 

320  FC0N=-1 

330  FOR  K=1  TO  3 

340  FOR  J=1  TO  3 

350  FC0N=FC0N+G<K>  J)*E*.N 

. J)*E<N,K> 

360  NEXT  J 

370  FCON*FCON+S<K>*E-'N,K 

> 

380  NEXT  K 

390  IF  FCON>0  THEN  G$=“F 
AIL":  RETURN 

400  IF  FC0N< -E5  THEN  440 

410  NC=NC+1 

420  CX':NC,1)=P 

430  CV.*NC,2)=N 

440  NEXT  N 

445  NEXT  P 

450  RETURN 


290  -  if  ply  thickness  zero, 
ignore  constraint 


300  -  Get  G  matrix  for  ply  being 
tested 


320  -  380  Solve  FCON  = 
Gijeiej  +  GiEi  "  1 


410  -  430  If  FCON  is  close  to 
zero  identify  constraint  as 
active,  make  list  in  C%  and  in¬ 
crement  constraint  counter 
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460  '**  GRAD IENV  ** 

475  UNORM=0 
480  II=P:  GOSUB  1230 
490  FOR  L=1  TO  NPLV 
508  IF  H(L>=0  THEN  700 
518  II=L:  GOSUB  1120 
520  FOR  J=1  TO  3 

539  Rt'.T>=0 

540  FOR  K=1  TO  3 

550  R< J>=R< J9-GKJ*  K>*E<N 
>  K  > 

560  NEXT  K,J 

570  FOR  J=1  TO  3 

530  Y< J>=0 

590  FOR  K=1  TO  3 

600  V< J>=VCnP+AKJ>K>*R< 

K) 

610  NEXT  K, J 

620  Z<L>=0 

630  FOR  J=1  TO  3 

640  FOR  K=1  TO  3 

650  2  <  L  >  =Z  < L  > +G  CJ , K  >  *  < V< 

J  >  *  E  <  N ,  K  > + E  C  N  •  J  >  *  Y  <  K  )  ) 

660  NEXT  K 

6  70  2  <1 > =2 < L  >  +S  <  J  > * V < J  > 
680  NEXT  J 

690  UN0RM=UN0RM+2 < L  >*2 <  L 


) 

700  NEXT  L 

710  UNGRM=SQR  <  UNORM  > 

720  FOR  L=1  TO  NPLV 
730  2<L>=Z<L) xIJNORH 
740  NEXT  L 
760  RETURN 
770  ’**  STRAINS  ** 

780  DIM  F<3»3> 

790  FOR  1=1  TO  3 
300  FOR  J=1  TO  3 
S10  F<  I  >  J>=A<  I  >  J >+[•'•  I  >  J) 

jfC 

820  NEXT  T.I 

830  DET#=F  (1,1  >*F(2,2:>*F 
( 3  ?  3  > + 2*F  < 1 » 2>  *F  <  2  >  3 )*F  < 
1 , 3)-F<2> 2>*F> 1 , 3>*F< 1 , 3 
)-F<l» 1>*F<2»3>*F<2,3>-F 


<3.3>+Fa,2)*F<l,2> 

340  Aia,l>=<F(2,2>»F(3, 
3>-F<2»3^*F<2j  3>  VOET# 
850  AI<2»2)=<F<1,1)*F<3> 
3)-F< 1 >  3)#F( 1 >3>>/C€T# 


360  AI<1,2)=<.F<1,3)*F<2, 
3>-F< 1 , 2)*F<  3. 3> )xDET# 
870  AI<3,3>=<F<1, 1 
2>-F<l,2>*F<1.2>VDET# 
830  AI<1»3)=<F(1>2>*F<2, 
3)-F<2>  2)*F< 1 >3>  >xDET# 
890  AI<2,3)=«:Fa,2>+F<l, 
3)-F< 1 > 1 >*F<2»3> >xDET# 


900  AI<2,l>=AI<l,2'oAI<3 
>2)=AI<2>3): AI<3> 1 >=AI<1 
.3) 

910  ERASE  F 


920  FOR  1=1  TO  NL 

930  FOR  J=1  TO  3 

940  E< I » J>=0 

950  FOR  K.=1  TO  3 

960  E<I.J>=E<I.J>+AIO.K 

>*XN<I,tO 

970  NEXT  K.J.I 

980  RETURN 


480 

ply 


-  Get  G  matrix  for  designated 


510  -  For  each  ply,  get  Q  matrix 
540-560  R  -  -  ~  A  e 


580  -  610  Y  =  A 


-1 


y  =  —1 _  e 

3hk 


620  -  680  $  (FCON)  =  [G..(e. 


ij'  i  3h, 


9  e  i 

3h, 


Ej»  -  Gi  ‘Ikt"  ik 


‘k  J  *  ““k 

690  -  730  Normalize  V(FCON) 


790  -  820  "F"  is  the  A  matrix 
corresponding  to  a  point  S 
along  the  Z  vector 

830  -  900  invert  A 

920  -  S70  Solve  e=|A  |  N  for 
each  independent  loading 
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990  ’**  A  MATRIX  ** 

1000  FOR  1=  1  TO  3 
1010  FOR  J=1  TO  3 
1020  AC  I , J)=0!  D< I >  J)=0 
1030  NEXT  J» I 
1040  FOR  1=1  TO  NPLY 
1050  I I=ls  gosub  1120 
1060  FOR  J=1  TO  3 
1070  FOR  K=1  TO  3 
3080  A>:jiK')=A<J»K)+Q<Jf K 
)*h<:  i ) 

1090  D< J»K>=D<J>K>+Q<J.K 
)*Z< I ) 

1100  NEXT  K»  J» I 
1110  RETURN 
1120  ’**  FORM  Q  ** 

1130  GK1, 1)=C(II. 1> 

1140  Qa,2)=C<I!»3> 

1150  QC1..3)=C(II,5) 

1170  Q<3,1>=C<II,5> 

1130  Q<3,2>=Cai,6> 

1190  GK3.3)=C(II»4) 

1195  CK2>3)=C<II»6) 

1208  GK2»2>=C<II»2) 

1210  GK2> 1 >=C<II>3) 

1220  RETURN 
1230  ’**  FORM  G  ** 

1248  6(1,1)-B(II>1) 

1250  6<lf2)*Bai»3) 

1260  sa»3>=e<:ii*5> 

1270  G<2,  l>=BdI»3> 

1280  GC2.2>=B<II»2> 

1290  G<2j3)=B<II j6> 

1300  Gv3- 1 >=B< 1 1 ? 5) 

1310  G<3,2>=B<II,6> 

1320  G<3,3)=B<II»4) 

1330  sa>*eai,7> 

1340  S<2)=B(II,'3> 

1350  S<3)=B<II,9> 

1360  RETURN 


1000  -  1100  The  matrix  D  is 
formed  so  that  along  the  Z  vector 

|A|  =  | A I  +  | D I  •  S 

where  S  is  a  scalar 


1130  -  1210  Convert  C  array  into 
3  x  3  Q  matrix  for  ply  designated 
by  II 


1240  -  1350  Convert  B  array  into 
3  x  3  G  matrix  for  ply  designated 
by  II.  Linear  failure  terms 
placed  in  vector  S 


1370  ’»*  HEW  H  UECTQR  ** 
1380  SMAX=1E10 
1390  FOR  1=1  TO  NPLV 
1400  IF  Z<I><>0  THEN  S=- 
H< I >/Z< I ) 

1418  IF  S>0  AND  S< SMAX  T 
HEN  SMAX=S 
1420  NEXT  I 
1438  Ft*"" 

1448  IF  SMAX>  10  THEN  Ft 
="FAIL"!  RETURN 
1450  S1=0!  S2=SMAX:  S=SM 
AX 

1460  IF  NC=0  THEN  1590 
1470  GOSUB  770!  GOSUB  26 
0 

1480  IF  G#="FAIL"  THEN  S 

2=S  ELSE  S1*S 

1498  IF  S1=SMAX  THEN  153 

1500  S=<Sl+S2)/2 
1510  IF  S2-S1<E2  AND  SI* 
0  THEN  F$="FAIL"!  S=0!  G 
0T0  1650 

1520  IF  Sl/'<S2-Sl)<4  THE 
N  1470 


1380  -  1420  Find  distance  along 

Z  to  find  h.  =0  constraint 
x 


1450  -  1500  Bisection  method  to 
find  distance  to  next  constraint. 
If  no  constraints  violated  at 
S  =  SMAX  then  stop  search 
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1530  S=S/2 
1535  SREF=0 
1540  FOR  1=1  TO  NPLV 
1550  H< I>=H< I >+Z< I >*S 
1560  IF  H<IXE2  THEN  H(I 
>=0 

1570  SREF=SREF+H< I >*H< I ) 
1530  NEXT  I 

1530  S=0 s  SREF =SQR  <  SREF  > 

160O  GOSUB  330:  GOSUB  77 

0:  GOSUB  2020 

1610  IF  SREF-S<E2  THEN  F 

*=BFAIL" 

1620  FOR  1=1  TO  NPLV 
1630  H<.  I  >=H<  I  >*S/SREF 
1640  NEXT  I 
1650  S=0 

1660  GOSUB  330:  GOSUB  77 

0s  GOSUB  260 

1670  RETURN 

1630  ’**  DIRECTION  ** 

1690  2=0:  IJNORM=l 

1700  FOR  1=1  TO  NPLV 

1710  X<I)=0 

1720  2=2+SGN<H(I)) 

1730  NEXT  I 
1740  Z=1/SQR<Z> 

1750  IF  NC=0  THEN  1360 
1760  FOR  1=1  TO  NC 
1 770  P=C:-.< I.  1 ):  N=CX< 1 . 2 
) 

1730  GOSUB  460 

1790  FOR  J=1  TO  NPLV 

1800  LET  X< J>=X< J>-Z<  J> 

1310  HEX.  M 

1815  UHORM-O 

1820  FOR  J=1  TO  NPLV 

1 830  UNOPM=UNCRM+X  < J ) *X< 

J> 

1340  NEXT  J 

1850  UN0RM=3QR ( UNORM ) :  T 
EST=0 

I860  FOR  1=1  TO  NPLV 
1370  X< I)=X< I ) /UNORM 
1830  TEST=TEST+X< I >»Z*SG 
N<H< I ) ) 

1890  NEXT  I 

1900  UNORM=0 

1910  FOR  1=1  TO  NPLV 

1920  Z<I)=X<I)-TEST*Z*SG 

N<H(I» 

1 930  UNORM=UNORM+Z < I ) *Z ( 
I) 

1940  NEXT  I 

1950  IF  UNORMC1E-6  THEN 
F*="F0IL":  RETURN  ELSE 
F$="" 

1960  UNORM=SGR ( UNORM ) 
1970  FOR  1=1  TO  NPLV 
1930  Z<I>=Z<IVUNORM 
1990  NEXT  I 
2000  GOSUB  990 
2010  RETURN 


1530  -  1600  at  point  halfway 
between  constraints,  use  strain 
ratio  routine  to  find  how  much 
the  laminate  thickness  can  be 
reduced 


1610  If  change  in  thickness 
small,  set  flag  to  halt  program 


1620  -  1660  Update  h  vector,  A 
matrix,  strains 


1760  -  1840  For  each  active 
constraint  call  gradient  sub¬ 
routine.  Sum  negative  of  each^ 
gradient  into  X  and  normalize  X 


1860  -  1890  Take  dot  product  of 
X  and  unit  normal  to  Ehi  =  const, 
plane 


1910  -  1940  Z  is  a  vector  parallel 
to  the  Zhi  =  const,  plane  and 
pointing  away  from  the  active 
constraints 


1950  if  the  magnitude  of  Z 
is  very  small,  a  local  minima  has 
been  reached 
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2028  '**  STRAIN  RATIO  ** 
2830  FOR  P=1  TO  NPLV 
2040  IF  H<P)=0  THEN  2160 
2050  II=P:  GOSUB  1230 
2060  FOR  N=1  TO  NL 
2070  B#=0*.C#=0 
2080  FOR  1=1  TO  3 
2090  FOR  J=1  TO  3 
2100  C#=C#-SREF*SREF*G<I 
, J>*E<N, I)*E(N, J> 

2110  NEST  J 

2120  B#=B#-SREF*S<I>*E<N 
.  I) 

2130  NEST  I 

2140  SUAL=<-B#+SQR  <B#*B 
#-4*C#  *  <1 -E6  >  >  > / <  2*  <1 -E6 
>) 

2150  IF  SUAL >S  THEN  S=SU 
AL 

2155  NEST  N 
2160  NEST  P 
2180  RETURN 
2190  »**  IFP  ** 

2200  2=1-SQR<NPLV) 

2210  FOR  1=1  TO  NPLV 

2220  2CI)=Z:  H(I)=Z 

2230  NEST  I 

2240  GOSUB  990 

2250  S=0:  SREF=1 

2260  GOSUB  770!  GOSUB  28 

20 

2270  FOR  1=1  TO  NPLV 
2280  H<n=H<I)*S 
2290  NEST  I 
2300  S=8 

2310  GOSUB  990!  GOSUB  77 

0:  GOSUB  260 

2320  RETURN 

2330  ’**  TRANSFORM  •** 

2340  FOR  1=1  TO  NPLV 

2350  C2=C0$  (.  2*T <!)’»!  C4= 

COS<4*TO>> 

2360  S2=SIN<2*T<I)>!  S4= 
SIN<4*T<I>> 

2370  B< I > 1 )=U< 1 )+C2*U(2) 
+C4*U<3> 

2380  B(I,2)=Uf 1)-C2*U<2) 
+C4*U(3> 

2390  B< I j3>=U<4)-C4*U<3) 
2400  e<I,4>=U<5>-C4*V<3> 
2410  B<  I  »5)=S2^2*U<'2)+S4 
*U<3> 

2420  B <  I > 6 )  =S2-'2*U <2'> -S4 
*U(3) 

2430  B< I »7)=U<6)+C2*U<7) 
2440  e<I,3>=<vK6)-C2*U<?> 
2450  B( I -9>=S2*U<7> 

2460  C< I » 1 )=U( 1 >+C2*U<2> 
♦C4*U<3> 

2470  C< I >2)=U<1)-C2*U<2> 
+C4*U(3> 

2480  C<I,3)=U<4>-C4*U<3> 
2490  Ca»4)=U<5)-C4*IJ-:3> 
2500  C<I,3)=S2/2*U''2>+S4 
*U<3) 

2510  C<I.6)=S2/2*U<2)-S4 
*U<3> 

2520  NEST  I 
2530  RETURN 


2030  -  2140  For  each  possible 
constraint  solve  for  S  in 


(SREF) 
e,E,  - -p- 


1]  i  ] 


+  G .  e  . 


l  l 


(SREF) 

S 


-  1  =  -E6 


2150  Take  smallest  value 
(corresponds  to  closest  constraint) 


2200  -  2310  For  equal  ply  ratios, 
find  the  smallest  laminate  thick¬ 
ness  which  does  not  violate  any 
constraints.  Initialize  A  matrix, 
strains,  and  constraint  list 


2370  -  2450 

Transform  failure 

parameters  in 

following  order 

B(I,1)=G11 

B(I,5)=G16 

B(I,2)=G22 

B(I,6)=G26 

B(I,3)=G12 

B(I,7)=G1 

B(I,4]=G66 

B(I,8)=G2 

B(I,9)=G3 

2460  -  2510 

Transform  modulus 

in  following  order 

C(I,1)  =  Qn 

C ( I , 5)  =  Q16 

c(i,l)-q22 

C(I,6)=Q26 

C(I,3)=Q12 

C ( I , 4 ) =Q66 
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2540  ’***  INPUT  *** 

2550  CLS 

2600  PRINT  "PRESS  ANV  KE 

V  WHEN": PR I NT "DESIRED  MA 
TEPIAL ":PRINT"AFPEhRS" 
2610  FOR  K=1  TO  750:  NEX 
T 

2620  FOR  M=0  TO  6 
2640  IF  M=6  THEN  Mf="HE0) 
MATERIAL"  ELSE  GETXM.EX 
. EV , OX • ES>  TPLV, XT , VT . XC , 
VC  >  SS . Mi 

2650  CLS SPRINT  M$:S0UHD2 
0, 1 

2660  FOR  J=1  TO  200 
2670  IF  INKEV$<>" "  THEN 
27O0 

2675  NEXT  J,M 

2680  GOTO  2620 

27O0  IF  11=6  THEN  GOSUB  9 

000: GOTO  2680 

2705  CLS: PRINT  "g  "JM^;" 

2710  PRINT  "HOW  MANV 
2720  INPUT  " PLV  GROUPS"; 
NPLV 

2770  CLS:  PRINT"ENTER  PL 

V  GPOUP" 

2740  PR I NT " CP I ENT AT I OHS " 
2750  FOP  1=1  TO  200 
2760  next  I 
2770  CLS 

2780  FOR  1=1  TO  NPLV 
2700  PRINT  " PLV  " ;  I 
2OO0  INPUT  TCI) 

2S10  TCI  -RIPADCTCI)) 
2S20  NEXT  i 

2830  PRINT"ENTER  NUMBER 
OF" 

2S40  PRINT  " INDEPENDENT 
LOAD" 

2350  INFUT  "CONDITIONS"? 
ML 

FOR  1=1  TO  NL 
2910  CLS:  PRINT  "Load 
I?"  in  MPa." 

2920  I NPUT  " N 1  =  " ; XN Cl , 1 ) 
2930  INPUT  "N2=";XN< 1,2) 
2940  INPUT  ,,N6=";XH<:i.3) 
2950  FOP  J=i  TO  3 
2960  XN1'  I  •  J)=XN<  I .  J>*1E6 
2970  NEXT  J, I 
2980  RETURN 


2600  -  2675  List  available 
materials.,  Get%  is  an  HX-20 
command  to  get  data  from  a  non¬ 
volatile  RAM  file 
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2990  ’**  INUARIENTS  ** 

3050  W=1/'<1-UX*UX*EV/EX 

3060  QXX=uv#EX*lE9;  qvv=  3050  -  3280  Calculate  invariants 

uv*ev*iE9  for  use  in  transformations.  Note 

3070  qxv=uv*ux*EV* l E9 :  Q  that  some  variables  like  EX  and 

3080 Vk? > = ( 3+QXX+3 ♦qvv+2  EY  9et  reused,  so  their  value  may 

*QXV+4*QS)"'8  not  be  what  you  might  expect  after 

3090  ik 2 > = < QXX-Q W > /2  routine  is  called 

3 1 00  l J  <  3  > = <  QXX+Q W -2+-QXV 
-4*0SV8 

3110  IK  4  > = <  QXX+Q VV +6*QXV 
-4*Q$)  8 

3 120  LK  5 ) = ( QXX+Q VV -2*QXV 
+4+QS)-S 

3130  EX=1E-12'<XT*XC>:  £ 

V= 1 E- 1 2-' ( VT* VC)  s  ES=1E-1 
2/'<S3+SS) 

3140  FX=  C 1  /XT- 1  .'XC >  ✓  1  £6 : 

FV=<1/VT-1/VCV1E6 

3150  e::v=-sqr  <ex*ew2 

3 1 60  GXX=EX*6XX*QXX+2*EX 
V*QXX*QXV+EV*QXV*QXV 
3 1 70  GW=EX*QXV*QXV+2*EX 
V*QXV*QVV+E  V*QVV*QVV 
3 1 m  GXV=EX*QXX*QXV+EXV* 
t.  Q  XX * Q V V + Q XV  *QX V  > + E V *QX V 
♦  QVV 

3190  GSS=ES*QS*QS 
3200  GX=FX*QXX+FV*QXV 
3210  GV=FX*QXV+FV*QVV 
3220  U ' 1 >  =  < 3+GXX+3+GVV+2 
*GXV+4+GSS)/8 
3230  U ' 2)= < GXX-GVV > /2 
3240  U C 3) = < 6XX+GYV-2#6XY 
-4+GSS )  "3 

3258  V ( 4 > =KGXX+GVV+6*GXV 
-4*GSSV8 

3260  U <5)*< GXX+GYV-2+GXV 
+4*GSS><8 

3270  U<6)*<GX+GV>/2 
3230  ',K7>=<6X-GY>/2 
3290  RETURN 
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3380  ’**  OUTPUT  ** 

3302  SOUND  15,2sSOUND50. 
2 

3305  K$="Hit  any  key  to 
cont.  " :  U$=  "  MN.-'in  " 

3310  CLS:  TEST=0 
3320  FOR  1=1  TO  NPLV 
3330  TEST=TEST+H<I>:  HEX 
T  I 

3350  PRINT  "TOTAL  THICKN 
ESS=" 

3360  PRINT  TEST! "  m.  " 
3370  PRINT  USING  "###«.# 
#  Plies" ’ TEST/TPLV 
3375  PRINT  K*? 

3380  A#=INKEV$: IF  A*<>"" 
THEN  3380 

3390  IF  INKEV$=" "  THEN  3 
390 

34O0  CLS:PRINT"Press  V  i 
f  pr intout" » "of  displace 
d  r  esu 1  t " , " i s  des i r  ed.  P 
ress  N","if  not"; 

3410  FOR  1=  1  TO  8008 NEX 
T  I 

3415  Af=INKEVi: IF  A*<>" " 
THEN  3415 

3420  CLS:  RESTORE  6120 
3425  J=Q: Af= INkEVt 
3430  FOR  1=1  TO  8 
3440  READ  AT: CLS: PRINT : p 
PINT  AT: SOUND20, 1 


3445 

AT=INKEVT 

: IFAT=""  T 

HEN 

3445 

3450 

PRINT  AT: 

: FOR  KK=1 

TO  75: NEXT  KK 

3455 

IF  AT=" V" 

THEN  J=J+ 

1 :  C; 

j  -  i.,=i 

3460 

NEXT  I 

3464  IF  JOG  THEN  LPRINT 
STRINGTC24,  "if'> 

3465  FOR  k:=  1  TO  J 

3470  ON  CM<K,1>  QOSUB  50 
00, 5200, 4000,4200,  4400, 4 
600 • 4800. 8O0O 
3485  LPRINT 
3490  NEXT  K 

3495  CLS: PRINT "FINISHED" 

,  KT 

3496  IF  INK.EVf =" "THEN349 
6  ELSE  RUN 
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4000  ’**  PLV  RATIO** 

4062  CLS: LPRINT  “Total  t 

hickness=" 

4004  LPRINT  US IMG  ".#### 
m.  "  >  TEST 

400c.  LPRINT  USING  "####. 
##  PI ies" ’ TEST/TPLV 
4O03  LPRINT 

4030  Af=" ANGLE  RATIO  « 
FI  IES" 

404O  LOCATE  0,1:PRINT  A* 
: LPRINT  AT 

405O  FOR  1=1  TO  HPLV 
4060  A=C I NT <  < FNDE6 < T< I > > 
+  1E2>'>  '1E2 

4m70  B=CINT < <H<I  VTEST*1 
E4  1  •  - 1E4 

4030  C=CINT (  <H(  I  VTPLV*1 
E 2>>'  1E2 

4030  PRINT  A;TAB<6)JB;TA 

b  <  1 3  ;  c 

4100  LPPINT  A;TABC6);B;T 

AB  C  13)?  C 
4120  NEXT  I 
4150  RETURN 

420O  ’**  STRENGTH** 

4210  LPRINT  "STRENGTH  RA 
TIOS" 

4215  LPRINT  "1=ULTIMATE 
STRAIN”: 

4220  LPRINT  ">1  IS  SAFE" 
4225  FOP  1=1  TO  NL 
4230  LPRINT  “LOADING  "I 
4235  LPRINT  "PLV" 

4240  FOR  P=1  TO  NPLV 
4245  IF  H<P)=0  THEN  4305 
4250  I I=p: GOSUB  1230 
4255  h#=0: E-#=0 
4260  FOR  J=1  TO  3 
4265  FOR  K=1  TO  3 
4270  A#=A#+G( J»K)*E< I » J) 
*E  ■'  I  ?  K  > 

4275  NEXT  K 

4230  B#«B#+S<J)*E<I,J> 
4285  NEXT  J 

4280  A#=(-B#+SQRCB#*B#+4 

*  A# >  >.'C2*A#> 

4295  A=FIX<A#*1E4+.  5)- IE 
4 

4300  LPRINT  FNDEGCTCPOS 

TAB' 10)' A 

4305  NEXT  P,  I 

4310  RETURN 

440O  •  *+STPA INS** 

4410  LPRINT  TAB<4) * "LAMI 
NATE  STRAINS" 

4420  FOR  N=t  TO  NL 
4470  LPRINT  "LOADING  "N 
4440  LPRINT  USING  "el=+# 

. ###E-03";E<N> 1)*1E3 
44501  LF'PINT  USING  "e2=+# 

.  ###E-03”;E'N,2)*IE3 
4460  LPPINT  USING  "e6=+# 

. ###E-03"?E<N,3)*IE3 

4465  NEXT  N 

4470  RETURN 

4600  ’  **A  MATRIX** 

4610  CLS 

4620  LPRINT "Norm.  |A|  in 
GPa.  " 

4630  FOR  1=1  TO  3 
4640  FOR  J=1  TO  3 
4650  DC  I. J)=A<I, J>- 1E9/T 
EST 

4660  NEXT  J,  I 
4670  GOSUB  7000 
4680  RETURN 


4210  -  4305  Strength  ratio  is 
defined  as  the  value  of  R  in 


G.  .  e  .  e  .R 
ID  i  D 


+  G. e .R  -1  =0 

l  l 


■  41.*  2  **  -.  •  '  iff  .'.V 
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4860  ’A  INVERSE 
4810  L PR I NT "Compliance 
< normal ized>" 

4820  LPRINT"  in  l  'TPa.  " 
4830  FOR  1=1  TO  3 
4840  FOR  J=1  TO  3 
4850  D(I,J>=AI<I,J?*TEST 

♦  1E12 

4 sea  next  j,i 

4870  GOSUB  7000 
4830  RETURN 

5000  LPRINT  "Material  Pr 
cperties" 

5O10  GETL'M- EX,EV,VX*ES,T 
PLV. XT i VT » XC-  VC*  SS? M# 
5015  LPRINT  M* 

5020  LPRINT  "EX=";EX?"GP 

5030  LPRINT  •,EV="5EV?"GR 
a" 

5040  LPRINT  "ES=";ES; “GP 
a" 

5050  LPRINT  "VX=";UX 
506O  LPRINT  "X="?XT; "MPa 

II 

5070  LPRINT  "X’="SXC?"MP 
a" 

5072  LPRINT  "Y="SVT;"MPa 

II 

5074  LPRINT  "V’s";VC;hMP 
a" 

5080  LPRINT  "S="?SS5 "MPa 

M 

5O90  LPRINT  "Pla  Thickne 
ss" !  TF’LV"m" 

5095  RE”'IRH 

600O  DA  In  10. 5E-5. .  1 . 1E- 
6 

6120  DATA  PI"  propertio 

*  Loads* Total  thickness 

ply  rati  os -Strength 
r at i os 

6130  DATA  Laminate  steal 
ns >  St i f i ness  m3tr i v » Comp 
liance  matrix* Ply  ratio 
araFh 

TOGO  ’ FANCV 

7010  LPRINT  - , — 

"  I  “  I  " 

702O  LPRINT  USING  "I###. 

1>*D<1*2>*D<1*3 

"l 

7030  A*="l - 1 - 

— t - P 

7040  LPRINT  AT 

7050  LPRINT  USING  " I ###. 

###"  ?  D  <2-1 .  ■ ,  D 2  -  2  >  -  D <  2 » 3 

> 

7060  LPRINT  A$ 

7070  LPRINT  USING  "I###. 
###"  5  D 3 » 1 ) , D<3> 2>  ■  l> <3-3 
) 

7O30  LPRINT  "  1 - 1 — 

_ i _ _ _ j  » 

7100  RETURN 
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8000  LPPINT  "**Ply  Ratio 
Gr  af=  h  *  ♦  " 

8010  LPPINT  "ANGLE" 

S02O  2=0: CLS 

802.0  FOR  1  =  1  TO  NPLV 

8040  n=H<o/TEST 
8050  if  ::<d;2  then  z=x< 


8040  :  I  -)=H<  I  >/TEST 
8050  IF  .  2  THEN  Z=X< 
I 

8055  NEXT 

8080  0ELTA=CINT<2--.  08>,  1 
00 

8070  A=CINT<DELTA*1000> 

808O  IF  AO20AND  A<  >50 AN 
D  h- ' 1O0AMD  a: >150  THEN 
DELI h=CFLTA+.  Ol: GOTO  807 
0 

8090  FOR  1=1  TO  NPLV 
8100  z='  X<I>.  DELTA*  12>+2 

f 

8110  I  I=CINT <ABS<COS<I -'2 
*3. 14159  > > > 

8115  IF  Z<26  THEN  81  GO 
8130  FOR  K=8  TO  15 
8140  H=t  *  1 1*  18: LINE C 28 , A 
> - • 2 »  A  >  ?  PSET 
8150  NEXT  K 

8  Ip.  O  LOCATE  0,1  + 1 1  *2:  PR  I 

NT  USING  "+###" 5 FN[>EG<T< 

I  :  LOChTEO,  0 

8170  IF  11=1  THEN  COPV:C 

L5 

81S0  NEXT 

8190  L PR  I NT  TAB<4) ? " I  I 
I  I  I  I  i  I  I  ” :  LF'RINT  TA 
B'6> ;  "1 " :  LPRINTU51NG" 

.  ##"?  DELTA:  LPRINT  TAE:< 
"Ft V  RATIO" 

8200  RETURN 

90O0  PRINT "REVIEW  OR  NEW 

M 

9010  INPUT  "RATA  ( FvN> " ; 
Af 

902O  IF  Ai="R"  THEN  9190 
9030  PR I NT" WHICH  MATERIA 
L":PRINT"UILL  VOIJ" 

9040  INPUT  "REPLACE  <0-5 
v :  i 

9050  INPUT  "EX<GPa>=";EX 
9O60  INPUT  "EV<GF'a  ; EV 
9070  INPUT  "UX=";UX 
9075  INPUT  "ESfGPa>=" ; ES 
9080  INPUT  "X(NPa)  =  " ;!< 
909O  INPUT  "X* <MPa>=" : XX 
9100  INPUT  “V<MPa)=": V 
9110  INPUT  " V ’ < MPa > = " 5 W 
9120  INPUT  "SCMPa)=";S 
9130  INPUT  "THICKNESS  <m 
)="? TPLV 

9140  INPUT  "NAME  <15  CHR 
.  MAX.  >  "  ?  Mi 

9 1 50  PUT7.I ,  EX ,  EV ,  UX >  ES ,  T 
PL V , X . V . XX , V V , S ,  Mi 
9160  PRINT  "ADDITIONAL": 
INPUT  “CHANGES  <Vx N>";Ai 
9170  IF  Ai="V"  THEN  9000 
9180  RETURN 

9190  PRINT"REVIEW  WHICH" 
:  INPUT " MATER  I AL  <0-5>" 5  M 
9200  GOCUB50G0 
9210  GOTO  9160 


1 


8020  -  8080  Automatic  scaling 


9150  -  Put%  is  an  HX-20  command 
to  place  a  data  string  into  a 
non-volatile  RAM  file 


vrv^-,M 
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9590  OPEN  MI",#1»  “CASl'-D 
ATto" 

9510  FOR  1=0  TO  5 
9523  INPUT  #1»EX, EV,UX,E 
S,  T,X,V,XX,VY,S,M* 

9530  PUT*: I ,  EX, EY,  UX ,ES,T 
,  X » V ,  XX ,  W«  S ,  M$ 

9540  NEXT 
9550  CLOSE  #1 
9560  DELETE  45 


9500  -  9550  Routine  to  read 
RAM  file  data  from  tape  called 
immediately  upon  loading  program 
and  has  no  further  function 


9560  -  Line  45  reads  GOTO  9500 
and  is  not  needed  after  data  is 
read 


5  DEFFIL  55,0 
10  ’**SAUE  RAM  FILE 
20  OPEN  "0",  #1,  "CASH DAT 
A" 

30  FOR  1=0  TO  5 

40  GET':  I ,  EX  -  EV,  UX>  ES,  T,  X 

, V, XX, VV' S  >  NT 

50  print  #i,ex;ev; nixes, 
t;xs v;  xxsvvs S?m$ 

60  NEXT 
70  CLOSE  #1 


9570  -  Routine  used  to  load  RAM 
file  onto  a  tape  after  program  has 
been  stored.  Usually  placed  in  a 
separate  program  area  from  main 
program  (using  command  LOGIN  2) . 
Only  needed  if  additional  tape 
copies  are  being  made. 


